Production of qq Pairs in Proton— Nucleus Collisions 

at High Energies 

Yuri V. Kovchegov a * and Kirill Tuchin b ' c ^ 
a Department of Physics, The Ohio State University, Columbus, OH 43210 
b Department of Physics and Astronomy, 
Iowa State University, Ames, I A 50011 
c RIKEN BNL Research Center, Upton, NY 11973-5000 
(Dated: February 2, 2008) 



We calculate production of quark-antiquark pairs in high energy proton-nucleus colli- 
sions both in the quasi-classical approximation of McLerran-Venugopalan model and 
including quantum small- x evolution. The resulting production cross section is ex- 
plicitly expressed in terms of Glauber-Mueller multiple rescatterings in the classical 
case and in terms of dipole-nucleus scattering amplitude in the quantum evolution 

n 

case. We generalize the result of |1|] beyond the aligned jet configurations. We expand 
on the earlier results of Blaizot, Gelis and Venugopalan by deriving quark pro- 
duction cross section including quantum evolution corrections in rapidity intervals 
both between the quarks and the target and between the quarks and the projectile. 
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I. INTRODUCTION 

Heavy quark production in hadronic collisions in high energy QCD is one of the most 
interesting and difficult problems. It is characterized by two hard scales: heavy quark mass 
m and the saturation scale Q s . The threshold for the invariant mass of the quark q and 
antiquark q production is 2m. Therefore, if m is much larger than the confinement scale 
h-QCD, it guarantees that a non-perturbative long distance physics has little impact on the 
quark production [3| making perturbative calculations possible (for a review see 0). 

Unlike the quark mass, which is a property of the produced quantum state, the satura- 
tion scale Q 2 S characterizes the density of color charges in the wave function of each of the 
colliding hadrons 0, 3, EJ • It increases as a power of energy and a power of atomic weight 
A 10, III |- At high energies and especially in reactions with heavy nuclei it becomes signif- 
icantly larger than the confinement scale. It is the saturation scale which makes the strong 
coupling constant small, a s (Q s ) <C 1, insuring applicability of the perturbative approach 
to all high energy scattering problems [^]. For all processes involving heavy quarks with 
momentum transfer of the order of Q 2 S ~ m 2 large saturation scale implies breakdown of 
the collinear factorization approach. The factorization approach may be extended by allow- 



ng 



12, 



;he incoming partons to be off- mass-shell. This results in conjectured /^-factorization 
13, Q|. Although the phenomenological applications of the /c^-factorization approach 



seem to be numerically reasonable already at not very high energies [15| its theoretical sta- 
tus is not completely justified. Like collinear factorization it is based on the leading twist 
approximation. However, at sufficiently high energies, higher twist contributions propor- 
tional to (Q s /m) 2n become important in the kinematic region of small quark's transverse 
momentum, indicating a breakdown of factorization approaches. 

The fact that the saturation scale at high enough energies and for large nuclei is large, 
Q s ^> Aqcd, combined with the observation that the typical transverse momentum of 
particles produced in pA scattering is of the order of that saturation scale, leads to the 
conclusion that Q s sets the scale for the coupling constant, making it small. This allows 
one to perform calculations for, say, gluon production cross section in pA collisions using 
the small coupling approach [16| . The same line of reasoning can be applied to heavy quark 
production considered here: the saturation scale Q s is the important hard scale making the 
coupling weak even if the quark mass m was small. Having the quark mass m as another large 
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momentum scale in the problem only strengthens the case for applicability of perturbative 
approach. 

Resummation of leading higher twist corrections to all orders have been performed be- 
fore in the Color Glass Condensate (CGC) /saturation framework for other observables not 
involving heavy quarks. The problem of gluon production in pA collisions at high energies 
was solved by resuming the contributions which are enhanced by factors of a 2 s A 1 ^ ~ 1 and 
a s y ~ 1, where A is the atomic number of the nucleus, and y is the rapidity variable [16]. 
Surprisingly, after resuming all such contributions to the single inclusive gluon production 
one recovers the /^-factorization formula ^| first suggested for the high parton density 
systems in Indeed, for large transverse momenta of the produced gluons, kx Q s , 
after neglecting all higher twist (Q s /kx) n corrections, one recovers the usual leading twist 
A>r-factorization. It was quite amazing that /c^-factorization for gluon production survived 
after resumming all twist s lip . However, kx- factorization fails for the double inclusive gluon 
production cross section |17| |. as well as for the inclusive quark production Instead a 
more complicated factorization picture emerges. 

Indeed the fact that the produced gluon transverse momentum spectrum in pA collisions 
obtained in [16] still diverges proportional to ~ l/k T in the infrared introduces logarithmic 
dependence of total gluon multiplicity dN/dy (integrated over all transverse momenta) on 
the infrared cutoff, raising questions about the applicability of the perturbative approach 
for calculation of that observable. However, while it is likely that Q s sets the scale for the 
running coupling even in dN/dy, a formal analysis of the scale of the running coupling is 
beyond the scope of this paper and is left for future research. Similarly, if one is interested 
in obtaining total quark multiplicity dN q /dy from the results presented below, one should 
strictly speaking view them as derived for quark production in deep inelastic scattering 
(DIS) (where the photon's virtuality Q plays a role of the infrared cutoff keeping the physics 
perturbative), which may also be applicable to pA collisions. 

Our goal in this paper is to calculate production of quark-antiquark pairs in high energy 
proton-nucleus collisions and in DIS both in the quasi-classical approximation of McLerran- 
Venugopalan model [9j (summing powers of a^A 1 ^ 3 ) and including quantum small- x evolu- 
tion (summing powers of a s y). We generalize the result of Q] for the single inclusive quark 
production beyond the aligned jet configuration. We derive the double inclusive quark and 
antiquark production. We expand on the earlier results of Blaizot, Gelis and Venugopalan 



by deriving a cross section that includes quantum corrections in the rapidity intervals 
between the quarks and the target (powers of a s y) and between the quarks and the projectile 
(powers of a s (Y — y)). (Here Y is the total rapidity interval, and y is the rapidity of the 
produced qq pair, with being the rapidity of the target.) We generalize the approach of jsjj 
by taking into account valence quark rescatterings in the nucleus in the quasi-classical ap- 
proximation, and also by including the quantum evolution corrections. In the quasi-classical 
limit our result should be equivalent to solution of the Dirac equation in the background of 



classical fields, similar to the one performed numerically in |2l| for a collision of two nuclei. 

The paper is structured as follows. We will first derive the qq production cross section 
in the quasi-classical approximation in Sect. II. We will proceed by including quantum 
evolution corrections in the obtained cross section in Sect. III. We will conclude in sect. IV 
by discussing phenomenological applications of the obtained results. 



II. INCLUSIVE CROSS SECTION IN THE QUASI-CLASSICAL 

APPROXIMATION 

The diagrams contributing to quark-antiquark pair production in the quasi-classical ap- 
proximation are shown in Fig.^ The graphs shown in Fig. ^ are dominant in the light-cone 
gauge of the proton. The first diagram corresponds to incoming valence quark emitting a 
gluon, which splits into a quark-antiquark pair before the system hits the target. The second 
diagram corresponds to the case when the valence quark first emits a gluon, after which the 
system rescatters on the target nucleus, and later the gluon splits into a quark-antiquark 
pair. The third diagram corresponds to valence quark rescattering on the target nucleus, 
after which it produces a gluon which splits into a quark-antiquark pair. 

The calculation of the diagrams in Fig. ^ will proceed along the lines outlined in [22! I23I] 
(see [24] for a review) using light cone perturbation theory |25(. In coordinate space the 
diagram contributions factorize into a convolution of Glauber-Mueller multiple rescattering 
with the "wave function" parts, which include splittings q v — > q v G and G —>■ qq. 

We begin by calculating the "wave-function" parts. In each of the diagrams in Fig. ^ 
they correspond to the two-step splitting q v — > q v G — > q v qq. However, the fact that the 
splittings take place either in initial or final states depending on the diagram modifies the 
energy denominators, making the "wave-function" parts different in all three graphs. We 
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proton 




FIG. 1: Three main contributions to quark- antiquark production in the quasi-classical approxima- 
tion. 

will denote these "wave-function" parts ^2 and ^3 correspondingly, as shown in Fig. 
The calculation of ^2 and ^3 proceeds according to the rules of light cone perturbation 
theory (LCPT) |2£j in the light cone gauge of the proton, which we choose as moving in 
the light-cone "plus" direction (see Fig. QJ). Calculations are first performed in momentum 
space, after which the "wave-functions" are Fourier-transformed into coordinate space. 

The important subtlety of calculating final-state splittings is that the light cone denomi- 
nator for such splittings should be calculated subtracting the light cone energy (the "minus" 
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momentum component) of the outgoing final state. Indeed the light cone energies of incom- 
ing and outgoing states are equal to each other: therefore, in calculating final state splittings 
one can still subtract the incoming energy in the denominators. However, in doing so one 
has to keep track of a change in the minus component of the target's momentum, which 
could be a bit tedious. For details on calculations of final state emissions in the LCPT 
formalism see 17 1. 

Since eikonal multiple rescatterings do not change the transverse coordinates of the in- 
coming quarks and the gluon, we can calculate ty 2 and ^3 in transverse coordinate space 
by calculating the diagrams in Fig. without interactions. We assume that the outgoing 
quark and anti-quark have momenta k\ and k 2 correspondingly. The "plus" components of 
the momenta, k\ + and k 2+ are conserved in the interactions with the target. Therefore, the 
"plus" component of the gluon's momentum is equal to k± + + k 2+ . Here, for simplicity, we 
assume that k i+ ,k 2+ <C p + , where p + is the typical light cone momentum of the valence 
quarks in the proton. This implies that k\ + + k 2+ <C p+, i.e., that the gluon is also much 
softer than the proton. In this kinematics the "wave-functions" in momentum space are 

K!AkiM) = 2gT a ^ gT b -= 7 ■ e - 



A 



2 9 J-ahj. — r^— irrzi — rr— 7 n — -7+- 



o ™ i* A ■ (£1 + £2) ^ ^,^(^1 (!-«)- *a«;«) 



A 



(Ai + fc 2 ) 2 M(l-a) + fc|a + m 2 



- 4# 2 T T ^.^(i-op ^ 
k^ (1 — a) + k%a + m 2 ' 

tf(2) (j. h n _ 9nT V e* A -fe + fe 2 ) 1 ^(fci) A jVgj 
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vi>(3) (h t,\- 9 T ST - (- 1 + T (h+ + k 2+ ) uJM) x v*>(h) 

A fc 1+ fe 2+ fci++fc 2 + fel+ k 2+ V 1+ V 2+ 



2 TT 1 1 u a (h) v a \k 2 



+ 2 ^ TaTb h + + k 2+ %^ + %^! v^7 7 " 

* A - (fci + fca) a (1 - a) _ ^ (1 - a) -.^a; a) 



A 



fc? (1 - a) + kj a + m 2 " \k x (1 - a) - k 2 a] 2 + m 2 



where 



k x (1 — a) + k 2 ot + m z 



a i{k x (1— a)— k 2 a; a) = — §. X -[ki (1— a)— A; 2 a] (1—2 a+A a) 5 CT)0 ./ -= a m (1— A a) ^-o-'j 

V2 

(2) 

A = ±1 is the gluon's polarization (which also does not change under eikonal rescatterings) , 
a = ±1 and a' = ±1 are quark and anti-quark helicities correspondingly (see Fig. ^ a' 
is defined with respect to — k 2 ), m is the mass of the quarks, and the colors of the gluon 
immediately after emission (a) and just before splitting into qq pair (b) are kept different 
since the color of the gluon is likely to change in interaction (for ^ 2 ), due to which the color 
factors will be calculated separately. Gluon polarization vector for transverse gluons is given 
by e A = (0, 0, e A ), with e A = (1 + i A)/ y/2. The fraction of gluon's "plus" momentum carried 
by the quark is denoted by a = ki + /(k 1+ + k 2+ ). The gluon "propagators" in diagrams 
and \l/ 3 of Fig. [T] have instantaneous (longitudinal) parts 25] , which account for the second 
(additive) terms in Eqs. (|Ta|) and (jlcjl. 

Note that, as can be checked explicitly using (JTJ), 

^ + *») + k) = 0, (3) 

indicating, of course, that there can be no emission without interaction. 

One may worry that since the gluon in the second graph of Fig. ^interacts with the target, 

and, therefore the interaction will depend on the transverse coordinate of this gluon, instead 

(2) 

of calculating W a a ,(ki, k 2 ) as shown above in (jlbjl . one should separately calculate q v — > q v G 
and G — > qq transitions in momentum space, and then separately Fourier-transform each of 
the results into coordinate space. However, this is not necessary, since the gluon's transverse 
coordinate is uniquely fixed by the transverse coordinates of the quark x x and the anti-quark 
x 2 and by a (see e.g. (20, 27, 3]). The gluon's transverse coordinate is 

u = ax 1 + {1 — a) x 2 . (4) 
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If we perform the calculations for q v — > q v G and G — > qq splittings independently, and 
Fourier-transform each of them into coordinate space, the G — > qq component will come 
with a delta-function S 2 (u — ax 1 + (l — a)x 2 ), which vanishes after integration over u (which 
is an internal variable and has to be integrated over) fixing u at the value given by (@J). The 
result of this procedure is equivalent to a simple Fourier-transform of ^^ a i(ki, k 2 ) from (Jlbj) 
into coordinate space. 

The light cone "wave-functions" in transverse coordinate space are defined as 

d 2 k 1 d 2 k 2 



e-^-^^l+k^h), 2 = 1,2,3. 



(5) 



(2tt) 2 (2tt) 2 

Here we assume that the transverse coordinate of the valence quark which emits the gluon 
(which splits into a qq pair) is 0. 

To perform the Fourier transform of JHJ) it is convenient to introduce the following auxil- 
iary functions 

F 2 (x 1 ,x 2 ;a) = J dq Jx{qu) K x (x n \/m 2 + q 2 a(l - a)) \/m 2 + q 2 a(l — a) , (6) 
Fi(xnX 2 ;a) = J dqj^qu) Kq(x 12 ^ m 2 + q 2 a(l - a) J , (7) 
Fo{x x ,x 2 ;a) = J dqqj (qu) K ^xi 2 V 'm 2 + q 



'a(l — a 



(8) 



where U — 1 7/ 1 . 3^12 — 1 —'2 "> *^ 12 — —12 

I, and q = k t + k 2 . In terms of the functions defined 

in ©, © and © we obtain 



^Lix^x^a) 



2g 2 T a T b 
(2n) 2 



E 



F 2 {x 1 ,x 2 ;a 



—12 ' i 

X12 



[l-2a + Xa)5 a ^ 



+ Fi(x x , x 2 , a) —= a m (1 — \ a) 5 a - a r 
V2 



u ■ e 



*A 



U 



- 25 a>(T > a(l- a) F (x v x 2 ; a) >, 



(9a) 



,(2) , ^ _ 2 g 2 T a T b 

(2vr) 2 4- 



—12 ' i 

£12 



mKAmxyi) (1 — 2 a + A er) 5 aa i 



+ K (w X12) —j= am (1 — A a) 5 a> - a > 
v2 



m • e 



*A 



The last relation (}9~cj) follows from ©. 



(2) 



(9b) 
(9c) 
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Summation over A yields 




F 2 (x 1 ,x 2 ;a) 



x 12 u 



1 



\{l -2a)x 12 - u + ia Ui x 12j ] S a y 





where e i2 = 1 = —£21, e n = e 22 = 0, and, assuming summation over repeating indices, 
€ij Ui Vj = u x v y — u y v x . Also x\ 2 j denotes the jth component of the vector x 12 . 

Now that we have calculated the "wave-functions" in (jlUj) . we can proceed by calculating 
the qq production cross- sect ion. The relevant diagrams are shown in Fig. Eland are obtained 
by squaring the sum of contributions from Fig. ^ We will first calculate the parts of the 
diagrams in Fig. El which are due to the squares of the "wave-functions" from (jlOJ) . The 
resulting expressions will then be convoluted with the multiple rescattering parts of the 
diagrams. 

The qq radiation kernel is obtained by averaging the square of the sum of the "wave 
functions" from (fTU|) over the quantum numbers of the initial valence quark and summing 
over the quantum numbers of the final state quarks. Since we are interested, first of all, in 
the qq inclusive production cross section, where the transverse momenta of both the quark k x 
and anti-quark k 2 are fixed, in anticipation of a Fourier transform to transverse momentum 
space, we will keep the transverse coordinates of the quarks different in the amplitude and 
in the complex conjugate amplitude. Therefore, if the transverse coordinates of the quarks 
are x x and x 2 in the amplitude, we will denote them by y and y 2 in the complex conjugate 
amplitude, as shown in the first graph of Fig. El The result for the squares of the "wave- 
functions" is 



^ijix^x^y^y^a) = — ^a'i^v^ «) (ivH^ a ) > f > J = 1,2,3. (11) 



Here the sum over gluons' colors a and b simply implies a calculation of the color factors of 
the relevant diagrams, including traces over fermion loops. Indeed these color factors, while 
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FIG. 2: Diagrams contributing to quark-anti-quark pair production in the quasi-classical approx- 
imation. Disconnected i-channel gluon lines imply summation over all possible connections to the 
adjacent s-channel quark and gluon lines. 
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calculable in principle, are rather sophisticated, especially if we are interested in the double- 
inclusive qq production cross section. Therefore, to simplify the already quite involved 
calculations, we will calculate the diagrams in Fig. Elin the large- N c limit. The other reason 
for doing this is that, even though it is clear how to improve on the large- N c approximation in 
the classical limit, inclusion of quantum evolution beyond the large- N c approximation would 
involve the functional JIMWLK j3] evolution equation, a numerical solution of which is 
rather involved. Therefore, in the calculations of "wave-functions" squared below, the color 
factors will be calculated in the large- iV c limit. 

After a straightforward calculation we derive (we introduce the gluon's transverse coordi- 
nate in the complex conjugate amplitude v = a y +(l — a) y 2 with v = \v\ and y = y x — y 2 -, 

2/12 = kJ) 



112 

2 



^nix^x^y y ;a) =AC F (— j \ F 2 (x 1 , x 2 ; a) F 2 (y y ■ a) [(1-2 a) 

- 1 - 2 V n J I x i2 Vn u v 

x (£12 •!*) (y,o -v) + (tijUiX 12j ) {e k iv k y 12 i)) + F^x^ x 2 ; a) F x {y y ■ a) m 2 



+4 a 2 (1 - a) 2 F (x 1 ,x 2 ;a)F (y i ,y :2 -,a) - 2 a (1 - a) (1 - 2 



a 



x 12 u 



F2(x 1 ,x 2 ; a) 



y ■ y_ 1 ^ 

x F (y y ;a) + =^ — F 2 (y , y ; a) Fq^, x 2 ; a) \ , (12a) 
- 1 - 2 2/12 v - 1 - 2 J J 

^(x^x^y ,y ;a) = AC F ( — ) m 2 { K x {mx 12 ) K x (m y 12 ) [(1 - 2a) 2 

- 1 - 2 \tt J { x 12 y 12 u v 2 



11 ' V 1 

x {x 12 -u) (y u -v) + {eijUiXi 2j ) {e M v k y 12 i)] + K {mx 12 ) K {my 12 ) =^ > , (12b) 

$i2fe, x 2 ;y ,y ; a) = -4C F (—) m \ F^x^ x 2 ; a) K^my^) [(1 -2a) 2 

- 1 - 2 \ir J { x 12 2/12 u v 2 

U * V 

x {x 12 -u) (y -v) + {tijUiX l2j ) {e k iv k y l2 i)\ + m Fi{x x , x 2 ; a) K (my 12 ) =-= 

- 2 a (1 - a) (1 - 2 a) Fo(x 1; x 2 ; a) AT^m y 12 ) 1 , (12c) 

2/12 « 2 " J 
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^33^, x 2 ;y v y 2 , a) = x 2 ; y v y 2 ; a) + $ 22 (x x , x 2 ; y v y 2 , a) + ^{x^x^^y 



a] 



+ $2i(x 1 ,x 2 ;y v y 2 ;a) (12d) 
$13(^1,^2;^,^; a) = ~® u&i, V V V 2 1 ot) -® 12 (xi, £2; y v y 2 ;&) (12e) 
■^(^i,^;^,^") = -^zife,^;^,^;^) - $ 22 (x 1 ,x 2 ;y i ,y 2 ;a) (12f) 

$ij(£i,£2;y 1 >y 2 ; a; ) = ^(y^^i^a)- ( 12 s) 

Here Eqs. fll2dll2ell2t|) follow from (jlOcj) . (|12gD allows one to obtain $ 21 , $ 31 and $ 32 from 



(IT2ct . (H2i and pljl . 

Rescattering of g„, g„ G and g„gg configurations on a large nucleus brings in different 
factors, which we label Hjj depending on the diagram shown in Fig. El For the case of 
single-quark inclusive production cross section (when transverse momentum of one of the 
quarks is integrated over) they were calculated in As we mentioned above, the calcu- 
lations complicate tremendously for the double-inclusive qq production cross section we are 
interested in calculating here. We will, therefore, perform our calculations on the large-iV c 



limit. Introducing quark saturation scale 



23, 



24] 



A 2 

Ql = ^ % T(b) « 2na 2 s N cP T(b) (13) 

with p the nucleon number density in the nucleus and T(b) the nuclear profile function, we 
write 

Zxi(x 1 ,x 3 lV 1 >h> ( *) = e-^-^ a ^^Vl2i-iulAHfe^) a o;Mi/l2 3 -v a |A) ) (14a) 

S 22 (x 1 ,x 2 ;y 1 ,y 2 ;a) = e -|(^ 2 Q* in(i/|^| A) j (14b) 

Saafei, £2; i^,^; a) = 1, ( 14c ) 

Zni^x^y^y^a) = e~\ ( « )2 02 MVI*i-«l a)-| ( m ) 2 Q? Mi/|* 2 -,| A) ? (Md) 

SasCsx,®,;^,^; a) = e -^ 2 o?Mi/«A )) (14e) 

SisfeL.a.J^^a) = e -^?^MiMA)-i,lQ ?M i/, 2 A) (14f) 

with A some infrared cutoff. All other H^-'s can be found from the components listed in (|14|) 
using 

^iiix^x^y^y^a) = ^ji(y v y 2 ; x t , x 2 ; a) (15) 
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similar to ( |12gD . Note that in arriving at equations (|14|) we have used the fact that the 
valence quark rescatterings on the target cancel due to real-virtual cancellations in the first 
four graphs in Fig. El [22]. Such cancellations do not happen completely for a projectile 
dipole, as we will see in Section III. 

Using Eqs. (j!2)) and (jTj| we derive the double-inclusive quark-anti-quark production 
cross section in pA collisions in the quasi-classical approximation 

da 1 



d 2 k 1 d 2 k 2 dydad 2 b A(2ix) 



d 2 x x d 2 x 2 d 2 Vl d 2 y 2 e-^rfe-j/xM^-fe-^) 



3 

x ^v(^^y v y 2 ^ a ) -ijix^x^y^y^a). (16) 

Here y is the rapidity of the s-channel gluon, which splits into the qq pair. Since the quark 
and the anti-quark are most likely to be produced close to each other in rapidity, one can 
think of y as the rapidity of the quarks, b is the impact parameter of the proton with respect 
to the nucleus. 

The single inclusive quark production cross section is easily obtained from (|16|) by inte- 
grating over one of the quark's momenta 

3 



da 



^—p- j d 2 x l d 2 x 2 d 2 y 1 | rfae iki ^ t>) ^ i 3 {x^ x 2 ] y_ v 



d?kdyd?b 2(2 



xHijfe,^;^,^;^, (17) 

where we inserted an overall factor of 2 to account for both quarks and anti-quarks. In (|17|) 
y is the rapidity of the produced (anti-) quark. 



III. INCLUDING QUANTUM EVOLUTION 



Here we are going to include small- x nonlinear quantum evolution of |29| into the cross 
sections from Eqs. (fl"6|) and (|17jl . Since the evolution equations in [29] are written for the 
forward amplitude of a quark dipole on a nucleus, we have to first generalize (fTKJ) to the case 
of qq production in dipole-nucleus scattering. Indeed, strictly speaking our results would 
then only be applicable to particle production in deep inelastic scattering. However, our 
results below may still serve as a good approximation for gluon production in pA collisions 
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jl^l ]. The generalization of Eqs. ()16|) and (|T7j) to dipole-nucleus scattering is easily done 
by including emissions of the s-channel gluon in Fig. El by the quark and anti-quark in the 
incoming dipole. If the transverse coordinates of the quark and anti-quark in the incoming 
dipole are denoted by z and z x correspondingly with z 01 = z — z v we write 
da 



loi) 



J d\ l d 2 x 2 d 2 y 1 d 2 y 2 e- ik ^~ v -^ k ^- y -2 



d 2 k 1 d 2 k 2 dydad 2 b K ^ )U 4(2 



3 1 

x Yl Yl ( _1 ) fe+ ' $ yfe _ ^fe'^2-^;y 1 -^ ) y 2 -^;a) ^ij^x^^y^^^a), (18) 

i,j=l k,l=0 

where now we have 

Hufe,^,^;^,!^) = e -l^) a «MVk 1 -!£ 1 |AH(a-^) a 02MVl2 3 -« a |A) j (19a) 

SaafeL.^&J^Ka.^a) = fi^UMi/M^ (l 9b ) 

Sssfe,*^;^,^,^) = e -|4^MiA«A) ) (19c) 

Znix^x^z^y^y^z^a) = e~\ fe"^ 2 ^ ia(i/| a -«|A)-i (* 2 -*) 2 Q 2 in(i/|* 2 -,| A) ? (1Qd) 

S23fe,^^;^^ 2 ,a;«) = e-l^^^/I^IA), (19e) 
Hiafe,^,^;^,^,^;^) = e -| fe-a) 2 Q? Mi/I^-^IA)-! fe-^ ( ) 2 Q2 Mi/i^-^iA) _ ^ 

Again, all other Hy's can be found from the components listed in (JT§j) using 

(x^x^z^y^y^^a) = 'E. ji {y^y^,z l ',x 1 ,x 2 , z k ; a). (20) 



The inclusion of quantum corrections due to leading logarithmic fresummin g po wers of 
a s y) approximation in the large- N c limit is done along the lines of [n| (see also [ttI Isol 
and [24J for a review) using Mueller's dipole model formalism [32]. Since the integration over 
rapidity interval separating the quark and the anti-quark in the pair does not generate a 
factor of the total rapidity interval Y of the collision (i.e., does not give a leading logarithm 
of energy), the prescription for inclusion of quantum evolution is identical to the single gluon 
production case. We first define the quantity ni(z , z x \ w ;l£ij Y — y), which has the meaning 
of the number of dipoles with transverse coordinates Wq,w_i at rapidity y generated by the 
evolution from the original dipole z £u z 1 having rapidity Y. It obeys the dipole equivalent of 



the BFKL evolution equat 



ion |32 



33] 



dniiz^z^^w^y) a s N c f 2 z 01 

a z 2 



dy 2 7r 2 J 4o z 2i 



2 

ni{za,^m,yh.iV) + n iU 2 , z t ; w , w^, y) 
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(21) 



with the initial condition 



ni(z ,£ L ;w ,w 1 ;y = Q) = S(z - w ) - w^. 



(22) 



If the target nucleus has rapidity 0, the incoming dipole has rapidity Y, and the produced 
quarks have rapidity y, the inclusion of small- x evolution in the rapidity interval Y — y is 



accomplished by replacing the cross section from (|18|) by 
da 

d 2 ki d 2 k 2 dy da d 2 b 



16, 



24. 



(i er 



(23) 

Note that while the substitution in includes only linear evolution, it results from analyz- 
ing all the possible non-linear evolution corrections including all possible pomeron splittings 
between the projectile and the produced qq pair. As was originally shown in Q| the pomeron 
splittings cancel in the rapidity interval between y and Y, leaving only the linear evolution 
contribution included in (J23J). 

Inclusion of evolution in the interval between and y is accomplished by replacing the 
Mueller-Glauber rescattering exponents according to the following rule ^| 



e 4 



|fe --i) 2 Q^MVI-o-- 1 |A) ^ l-N{xn,Xz,Y) 



(24) 



where N(x ,x 1 , Y) is the forward amplitude for a quark dipole Xq,^ scattering on a target 
with rapidity interval Y between the dipole and the target. It obeys the following evolution 
equation |2Ql | 



dN^x^Y) _ a s N c 
dY 2 Tv 2 



d 2 x 2 -4^T l N (^ a, y) + N(x 2 , x^Y)- N^, x 1} Y) 

X 20 X 21 



w(^,x 2 ,r)iv(z 2 ,x 1 ,y)] 



(25) 



with the initial condition 



N(x ,x 1 ,Y = 0) = 1 - e -^o--i) 2 Q 2 s in(i/| £o -^i|A). 



(26) 
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Performing the substitution from ()24j) in (JT§j) yields 

^uix^^z^^y^z^^Y) = -Nix^y^.Y) - N(x 2 ,y 2 ,Y) 

+N{x 1 ,y 1 ,Y)N{x 2 ,y 2 ,Y), (27a) 

^22(^X2, z^y^y^^; a, Y) = -2 N(u,v,Y) + N(u,v,Y) 2 , (27b) 

^(x^x^z^y^y^z^^Y) = -2 N(z k ,^,Y) + N(z k ,z l7 Y) 2 , (27c) 
^i2{.Xi,x 2 ,z k ]y_ v y 2} z l ] a,Y) = -N(x x ,v,Y) - N{x 2 ,v,Y) 

+N(x 1 ,v,Y)N(x 2 ,v,Y), (27d) 

223(^1,^2,^; y^y^Zi'i^Y) = - 2N (^z u Y) + Niu.^Yf , (27e) 

2i3(x 1 ,x 2 ,z Jfc ;y 1 ,y 2 ,z i ;a,y) = _iV fe'^' r ) - ^(a^^r) 

+iV(x 1 ,^,^)iV(x 2) ^y) ) (27f) 



with 



Sijfe,x 2 ,2 fc ;y i ,y 2 ,^;a,r) = T E ii {y_ v y_ 2 ,z i \x x ,x 2 , z k ; a,Y). (28) 



In arriving at (j2"Tj) we have dropped additive unit terms which do not contribute to the 
cross-section due to (jHJ) leading to Y^Hj=i ®v = ®- 

With the definition of Eqs. (f2*Tj) we write the following answer for the double inclusive 
qq production cross section including small-x evolution effects 

%i) = ~rjn \r / d 2 WQd 2 w 1 n 1 {^z 1 ]w Jd ,w l -,Y -y) 



d 2 k l d 2 k 2 dydad 2 b y ^ JU 4(2vr) 6 

x d 2 Xl d 2 x 2 d 2 Vl d 2 y 2 e -^i-fe-Ii)-^2-fe-^ 2 ) 



3 1 

x Yl -^fc>^2 -Wjc',^ -m,y_ 2 ~m; a ) -iji^x^w^y^y^^a^y). 

i,j=l k,l=0 

(29) 

Similar to how we arrived at (|17|) . we integrate over one of the quarks' transverse momenta 
to obtain the single inclusive quark production cross section 

U01) = TTTo — \Z f d 2 w d 2 Wi ni(z Q} z^WqjUIi, Y — y) d 2 X\ d 2 X2 d 2 y\ e" 1 -^ 1 *-^ 



d?kdyd 2 b y -" u 2(2tt) 



»1 3 1 

^° i,j=l k,l=0 

(30) 



x 
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Equations and (jSOj) are the central results of this paper. 

An example of the pomeron fan diagram contributing to the obtained cross sections is 
shown in Fig. El There the proton or qq dipole in DIS is shown on top of the figure. The 
nucleus is represented by the straight lines at the bottom of the figure. As usual each ladder 
represents a BFKL pomeron. Fig. El reflects the main features of equations (J2H|) and (J3U|) : 
it contains a linear evolution between the produced qq pair and the projectile, and the non- 
linear evolution between the qq pair and the target. Indeed a relatively simple diagram in 
Fig. El does not include all the complexity of the nonlinear interactions in ()27j) and of the 
emission wave functions in fT2|) . 

i 
i 




FIG. 3: An example of a pomeron fan diagram contributing to the qq production cross section in 
proton- nucleus collisions or in DIS as calculated in equations and (|30j). The produced quark 
and anti-quark are marked by crosses. 

IV. SUMMARY 



Expressions ()29|) and (|30|) for the single and double inclusive quark production have been 
derived by summing perturbation series in the coupling constant a s . In that sense our 
result is perturbative. It was pointed out in [39|, |40|, |4l| that there can be a significant 
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non-perturbative contribution to particle production in high energy QCD. The investigation 
of this effect is beyond the scope of the present paper: however it certainly deserves further 
study 

Equations ()29j) and (jHOj) have important phenomenological applications for studying the 
dense partonic system in p(d)A and eA collisions. Observation of hadron suppression in the 
nuclear modification factor measured in dA collisions at forward rapidities at Relativistic 
Heavy Ion Collider (RHIC) |34| signals the onset of the nonlinear evolution of the scattering 



amplitude for light hadrons |35|, |36|, [j2|. Due to a large mass, the impact of nonlinear 
evolution effects on the heavy quark production is shifted to higher energy and/or rapidity. 
It was estimated in 19] using the /^-factorization approach that one can expect a significant 
deviation of the open charm production cross section from the perturbative behavior already 
at pseudo-rapidity r] ~ 2 at RHIC. Due to the heavy quark production threshold one expects 
that the total multiplicity of open charm scales as N co \\ at lower energy and/or rapidity 
whereas at higher energies and/or rapidities the scaling law should coincide with that for 
lighter hadrons [l^, i. e. open charm multiplicity should scale as iVp art due to high 
parton density effects. Therefore, to be able to compare predictions of CGC with the data 
reported by RHIC experiments and to make predictions for the possible upcoming pA run 
at the Large Hadron Collider (LHC), it is important to perform a calculation of an open 
charm production within the more general approach developed in this paper. Our final 
results (|29|) and (J30|) allow one to describe open charm transverse momentum spectra at 
different rapidities and center-of-mass energies, allowing for a complete description of RHIC 
and LHC data. Since the saturation scale Q s is expected to be even higher at LHC than 
it was at RHIC, the CGC effects on heavy quark production at LHC should be even more 
significant. 
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